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ABSTRACT 

We use the microlensing variability observed for nine gravitationally lensed quasars to show that 
the accretion disk size at 2500A is related to the black hole mass by log(i?2500 /cm) = (15.6 ± 0.2) -|- 
(0.54± 0.28) log{MBH /W^Mq). This scaling is consistent with the expectation from thin disk theory 

{R cx m"^^), but it implies that black holes radiate with relatively low efhciency, log(?7) = — 1.29± 
0.44 -I- \og{L/LE) where 77 = L/{Mc^). These sizes are also larger, by a factor of ^ 3, than the size 
needed to produce the observed 0.8 /xm quasar flux by thermal radiation from a thin disk with the 
same T cx R~^/'^ temperature profile. More sophisticated disk models are clearly required, particularly 
as our continuing observations improve the precision of the measurements and yield estimates of the 
scaling with wavelength and accretion rate. 

Subject headings: accretion, accretion disks — dark matter — gravitational lensing — quasars: general 
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I. INTRODUCTION 

Despite nearly 40 yea rs of work on a cc retion 
disk physics, the simple IShakura fc SunvaevI ()1973D 
thin disk model and its relativistic cous i ns (e.g. 
Page fc Thorndll974l : Hubenv fc Hube!^ [l997l : IlIetaLl 
2ooi r remain the standard model despite many theo- 
retical alternatives ( e.g. Naravan, Barret & McClintoc^ 



[19971 iDe Villiers. HawTev fc KrolikI 120031 : IBlacs 2007 



I Blacs 2007i) 

and some observational reservations (see ICollin et all 
|2002| ). Quasar accretion disks cannot be spatially 
resolved with ordinary telescopes, so we have been 
forced to test accretion physics through time variabil- 
ity (e.g. LV andcn Berk ct al. 2004'; 'Serg eev et all 120051 : 
ICackett. Hornc fc Winkler. .2007) and spectral m odeling 
fe.g. iSun fc MalkanI 119891: iBonning etall [200I . One 
notable success is t he use of reverberation mapping 
(jPeterson et al.ll2004D of quasar broad line emission to 
calibrate the relation between emission line widths and 
black hole masses. 

* Based on observations obtained with the Small and Moderate 
Aperture Research Telescope System (SMARTS) 1.3m, which is 
operated by the SMARTS Consortium, the Apache Point Obser- 
vatory 3.5m telescope, which is owned and operated by the As- 
trophysical Research Consortium, the WIYN Observatory which 
is owned and operated by the University of Wisconsin, Indiana 
University, Yale University and the National Optical Astronomy 
Observatories (NOAO), the 6.5m Magellan Baade telescope, which 
is a collaboration between the observatories of the Carnegie In- 
stitution of Washington (OCIW), University of Arizona, Harvard 
University, University of Michigan, and Massachusetts Institute of 
Technology, and observations made with the NASA/ESA Hubble 
Space Telescope for program HST-GO-9744 of the Space Telescope 
Science Institute, which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under NASA contract NAS 
5-26555. 

^ Department of Astronomy, The Ohio State University 



Gravitational telescopes do, however, provide the nec- 
essary resolution to study the structure of the quasar 
continuum source. Each gravitationally lensed quasar 
image is observed through a magnifying screen created 
by the stars in the lens galaxy. Sources that are smaller 
than the Einstein radius of the stars, typically ~ 10^^ 
cm, show time variable fluxes whose ampli tude is deter- 
mine d by the source size (see the review by I WambsganssI 
[200i ). Smaller sources have larger variability amplitudes 
than larger sources. In this investigation, we exploit 
the optical microlensing variability in nine gravitation- 
ally lensed quasar systems to measure the size of their 
accretion disks, and we find that disk sizes are strongly 
correlated with the masses of their central black holes. 

In we describe the monitoring data, the lens models 
we use based on Hubble Space Telescope (HST) images of 
each system and our microlensing analysis method. In fj3] 
we describe our accretion disk model and our results for 
the relationship between disk size and black hole mass. 
In 21 discuss the results and their implications for 
thin accretion disk theory. All calculations in this paper 
assume a flat ACDM cosmology with h = 0.7, Hm = 0.3 
and r^A = 0.7. 

2. DATA AND ANALYSIS 

We monitored the gravitationally lensed quasars 
HE 0435-1223, SDSS 0924-^0219, FBQ 0951-^2635, 
HE 1104-1805, PC 1115-F080, RXJ 1131-1231, SDSS 
1138-1-0314, SBS 1520-^530 and Q 2237-h030 in 
the R- and F-bands on the SMARTS 1.3m using 
the ANDICAM optical/infrared camera (jDePov et al.l 
l2003f )^. the Wisconsin- Yale-Indiana (WIYN) observatory 

^ http:/ /www. astronomy.ohio-state.edu/ANDICAM/ 
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TABLE 1 
Derived Quantities 



Object 


Mbh 


log(_Rs/cm) 




Icorr 


log(^?s/cm) 




(10^ Mq) 


(microlensing) 


(/im) 


(mag) 


(thin disk flux) 


HE0435-1223 


0.50 


15.6+0-5, 


0.260 


20.76 ±0.25 


14.9 ±0.1 


SDSS0924+0219 


0.11 


14.8t0i 


0.277 


21.24 ±0.25 


14.8 ±0.1 


FBQ0951+2635 


0.89 


15.9«;i 


0.313 


17.16 ±0.11 


15.6 ±0.1 


HE1104-1805 


2.37 


7+0.2 


0.211 


18.17 ±0.31 


15.4 ±0.1 


PG1115+080 


0.92 


16 5+° '5 


0.257 


19.52 ±0.27 


15.1 ± 0.1 


RXJ1131-1231 


0.06 


ir 4+0.3 


0.422 


20.73 ±0.11 


14.8 ±0.1 


SDSS1138+0314 


0.04 


14 8+°-'' 


0.203 


21.97 ±0.19 


14.6 ±0.1 


SBS1520+530 


0.88 


15 5+°-2 

10.0_Q 2 


0.245 


18.92 ±0.13 


15.3 ±0.1 


Q2237+030 


1.3 


ir 5+0.3 


0.208 


18.03 ±0.44 


15.5 ± 0.2 



Note. — Rs from microlensing is the accretion disk size at Arest, the rest-frame 
wavelength corresponding to the center of the monitoring filter used for that quasar's 
lightcurve. Icorr is the corrected (unmagnified) /-band magnitude. Typical /-band 
measurement errors are < 0.1 mag, but the larger errors on Icorr come from uncertainties 
in the lens magnification. Rg calculated using corrected /-band magnitude and thin disk 
theory is also unsealed; it is the disk size at the rest-frame wavelength corresponding 
to the center of the HST /-band filter (F814W). Both disk sizes assume an average 
inclination angle i = 60° . 



using the WIYN Tip-Tilt Module (WTTM) 3, the 2.4m 
telescope at the MDM Observatory usi ng the MDM 
Eight -K 4, Echelle and RETROCAM ^ (jMorgan et all 
I2OO5I) imagers and t he 6.5m Magellan Baade tele- 
scope using IMACS (jBigelow et alj |1999[ ). We sup- 
plemented our m onitoring data with p ublished quasar 
light curves fro m Paraficz et al. (2006*), Schechter et al.' 
([l997l).IWvrzvkowsk i e t al. (2003), Ofck & Maoz (2003) 
and iGavnullina et~an (|2005f) . We measured the flux of 
each image by comparison to the flux from reference stars 
in the field of each frame. Our analysis of the monito ring 
data is described in detail by iKochanek et al] (|200 6^. 

All nine lenses have been observed in the V- (F555W) , 
/- (F814W) and i7-bands (F160W) using the WFPC2, 
ACS/WFC and NICMOS instruments on HST. We fit 
these images as combinations of point sources for the 
quasars and (generally) de Vaucoul eurs models for the 
lenses as described in ()Lehar et al.l ,2000) . These pro- 
vided the astrometry used for lens models and defined a 
constant mass-to-light (M/L) ratio model for the mass 
distribution in the lens models. We model ed each system 
using the GRAVLENS software package ()Keetonl I2OOII) 
to generate a series of ten models for each lens starting 
from a constant M/L model an d then adding an NEW 

avarro. Erenk fc Whiti [l99l halo. The sequence is 
parametrized by /m/L; the mass fraction represented by 
the visible lens galaxy relative to a constant AI/L model. 
We start with the constant M/L model, /m/l ~ Ij and 
then reduce its mass in increments of /S.fM/L =0.1 with 
the NEW halo's mass rising to compensate. 

These lens models then provide the convergence k, 
shear 7 and stellar surface density needed to define the 
microlensing magnification patterns. We assume a lens 
galaxy stellar mass function dN{M)/dM oc M~^-^ with 
a dynamic range of a factor o f 50 that app roximates the 
Galactic disk mass function of lGouldl (|2000l ) . For the typ- 
ical lens we generated 4 magnification patterns for each 

^ http:/ /www. wiyn.org/wttm/WTTM_manual. html 

* http://www.astro.columbia.odu/ arlin/MDM8K/ 

5 http://www.astronomy.ohio-state.edu/MDM/RETROCAM 



image in each of the 10 lens models. We gave the mag- 
nification patterns an outer scale of 20{Re), where (Re) 
is the Einstein radius for the mean stellar mass (M). 
This outer dimension is large enough to fairly sample 
the magnification pattern while the pixel scale is small 
enough to resolve the accretion disk. We determined the 
properties of the accretion disk by modeling the observed 
light curves usin g the Bay esian Monte Carlo me thod of 
iKochanekl (|2004l ) (also see lKochanek et al.|[2007l ). For a 
given disk model we randomly generate light curves, fit 
them to the observations and then use Bayesian meth- 
ods to compute probability distributions for the disk size 
averaged over the lens models, the likely velocities of the 
observer, lens, source, stars, and mass. We use a prior 
on the mean microlens mass of O.IM© < (M) < I.OMq, 
but the disk size es timates are relat ively insensitive to 
this assumption (see 'Kocha nekll2003 ) . 

We use black hole mass estimates for the quasars 
that are based on observed quasar emission line widths 
and the loca lly calibrated viria l relations for blac k 
hole masses (|Onken et al.l [200I [G reene IHfol [200I . 
For most systems we simply ad o pted the black hole 
mass estimates from Pen g et all (|2006[ ) based on the 
C IV (AI549A), Mg II (A2798A) and (A486lA) mass- 
linewidth relations. For SDSS 1138+0314, we measured 
the width of the C IV (1549 A) hne in optical spectra from 
the S loan Digital Sky Survey (jAdelman-McCarthv et aH 
l2006( l and estimated the black hole mas s usin g 
the normalizations of IVestergaard fc PetersonI (j^OOG). 
These mass e stimates are reliable to approxi mately 
0.3 dex (see [McLure fc Jar vis 2002; KoUmeie r et al.l 
I2006t IVestergaard fc Petersonli200Q : iPeng et al.ii200l) . 

3. RESULTS 

We model the surface brightness profile of the accre- 
tion disk as a power law temperatur e profile, T oc R~^^^, 
matching the outer regions of a iShakura fc SunvaevI 
(1973) thin disk model. We neglect the central depres- 
sion of the temperature due to the inner edge of the disk 
and corrections from general relativity to avoid extra pa- 
rameters. The effect of this simplification on our size 
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Fig. 1. — Inclination-corrected accretion disk size -R2500 versus black hole mass Mbh- The solid line shows our best power-law fit to 
the data and the dot-dashed line shows the prediction from thin disk theory (L/Le = 1 and 77 = 0.1). Disk sizes are corrected to a rest 
wavelength of \rest = 2500A and the black hole masses were estimated using emission line widths. The filled points without error bars are 
^2500 estimates based on the observed, magnification-corrected 7-band fluxes. They have typical uncertainties of 0.1-0.2 dex. 



estimates is small compared to our measurement uncer- 
tainties provided the disk size we obtain is several times 
larger than the radius of the inner disk edge. We assume 
that the disk radiates as a black body, so the surface 
brightness at rest wavelength Xrest is 



2hr, 



\3 

^rest 



exp 



R 



3/4 



-I -1 



1 



(1) 



where the scale length 



9.7 X 10' 



hpC^ 



BH 



/im 

2/3 



1/3 



4/3 



/ L 



1/3 



loll , cm (2) 

is the radius at which the disk temperature matches 
the wavelength, kT\^^^^ — hpc/ Xrest, hp is the Planck 
constant, k is the Boltzmann constant, Mbh is the 
black hole mass, M is the mass accretion rate, L/Le 
is the luminosity in units of the Eddington luminosity, 
and 77 — L/{Mc^) is the accretion efficiency. We can 
also compute the size under the same model assump- 
tions based on the magnification-corrected /-band quasar 
fluxes measured in HST observations as 



ycos i V '"ff 

3/2 



cm 



(3) 



where Dos/ra is the angular diameter distance to 
the quasar in units of the Hubble radius, / is the 
magnification-corrected magnitude and i is the disk in- 
clination angle. 

Our results are shown in Figures [1] through [3] and 
summarized in Table [1] For the comparison with the- 
ory and the figures, we corrected the measured sizes to 
Xrest = 25OOA assuming the A'*/"^ scaling of thin disk 
theory and the mean inclination (cosi) — 1/2. There 
are two striking facts illustrated by the figures. First, 
we clearly see from Fig. [T] that the microlensing sizes are 
well correlated with the black hole mass. A power-law 
fit between i?2500 and Mbh yields: 



log ( (15.6 ±0.2) 



-(0.54 ±0.28) log 



/ M, 



BH 



VlO^Mc 



(4) 



which is consistent with the predicted slope from thin 
disk theory (i? cx m'^^) and implies a typical Eddington 
factor of log(L/?7L£;) = 1.29±0.44 if we fix the s lope with 
mass to 2/3 (see Fig. 2). iKollmeier et al.l (|2006[ ) estimate 
that the typical quasar has L/Le ~ 1/3, which would 
indicate a radiative efficiency of 77 = L/{Mc^) ~ 0.02. 
This efficiency is low compared to standard models (e.g. 
lGammidll999h . Second, we find that microlensing sizes 
are well correlated with sizes estimated from thin disk 
theory and from the observed flux (Figs. [T]and[2]), but 
the three size estimates show systematic offsets in scale. 
Most of the offset between the microlensing disk size 
measurements and thin disk theory size estimates could 
be explained by the existing uncertainties, but the off- 
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Fig. 2. — Results of the power-law fit to i?2500 as a func- 
tion of black hole mass. The contours show the 1 — 3cr one- 
parameter confidence intervals for the slope a and the normal- 
ization R250o(Mgu = IO^Mq) for the 2500A accretion disk size 
corresponding to Mgjj = IO^Mq. The best-fit value is indicated 
with a black point. The filled points along the dot-dashed line are 
theoretical thin disk sizes for quasars radiating at the Eddington 
limit and with efficiencies of ?7 = L/{Mc^) = 0.01, 0.1 or 1.0. 

set from the estimate based on the quasar flux is much 
more significant (Fig. [T|). While the microlensing and 
flux sizes are well-correlated, the measured disk sizes are 
0.4±0.3 dex larger than predicted from the observed flux. 
Simply put, the quasars are not sufhciently luminous to 
be radiating as black bodies with a T cx R'^^'^ tempera- 
ture profile. iPoolev et al.l (|2006f ) also noticed this prob- 
lem in their more qualitative study of lensed quasars with 
X-ray observations. 

4. DISCUSSION 

Are the discrepancies between our measurements and 
the size estimated from the disk flux due to a problem 
in the measurements, an oversimplification of the disk 
model or a fundamental problem in the thin disk model? 
We have tested our approach using Monte Carlo simula- 
tions of light curves and verified that we recover the input 
disk sizes. Our results are also only weakly sensitive to 
theassumed prior on the microlens masses (see Kochanel^ 
|2004 for a discussion). We will overestimate the source 
size if a significant fraction of the observed flux comes not 
from the continuum emission of the disk but from the 
larg er and minimally microlensed line emitting regions 
(e.g. lSugai et al.ll2007l ). This includes not only the obvi- 
ous broad lines but also the broad Fe II and Balmer con- 
tinuum emission that can represent ~ 3 0% of the appar- 
ent continuum flux at some wavelengths (|Netzer fc Willis! 
Il983t lGrandilll9 82fl. We have experimented with adding 
a fraction of unmicrolensed light and found that 30% 
contamination does not lead to sufficiently large size 
changes to resolve the problem. The sizes shrink by ap- 
proximately 20%. Conservatively, black hole mass esti- 
mates fro m the viri al technique have a scatter of a fac- 
tor of ^ 3 (iMcLure fc Jarvis 2002; KoUmcicr et al. 2006; 
IVestergaard fc Peterson! 120061 : IPeng et all 120061 ) . which 
contributes only 0.3 dex of scatter to the disk size es- 
timates. The size estimates from the flux could be af- 



FlG. 3. — Thin disk flux size estimates versus accretion disk 
sizes from microlensing. For reference, the solid line indicates a 
one-to-one relationship between thin disk flux size estimates and 
microlensing measurements. The dot-dashed line is the best fit to 
the data. Since the data points have large errors relative to their 
dynamic range, the best-fit slope is consistent with unity and its 
average offset from the solid line is 0.4 dex. 

fected by misestimating the magniflcation or failing to 
correct for extinction in the lens galaxy, but the uncer- 
tainties in the magnifications are only a factor of ~ 2 at 
wors t and none of these lenses shows signific ant extinc- 
tion (|Falco et al.lll999tlElfasd6ttir et al.ll2006l ) relative to 
the magnitude of the discrepancy. 

The problem is also not a consequence of the obvious 
flaws in our simplified disk model - the neglect of the 
inner edge and t he different temperatu re profile of a rela- 
tivistic disk fe.g. lPage fc Thornelll974l ). For observations 
at a fixed wavelength, neglecting the inner edge does not 
have a dramatic effect on the effective source size and 
the effects of relativity on the temperature profile are 
modest. Studies of the temperature profile in the disk 
with microlensing are best done by measuring the vari- 
ation in source size with wavelength because size ratios 
can b e measured much more accurately than absolute 
sizes. iPoindexter et al.] ()2007f ) have used the wavelength 
dependence of the microlensing in HE I104-I805 to de- 
rive a slope T cx R-/^ or Rx cx A^/^ of /3 = OMtoAl that 
is consistent with thin disk theory but would also allow 
a shallower temperature profile that would reduce the 
differences between the microlensing and flux size esti- 
mates. The next step is clearly to use more sophisticated 
disk mo dels including relativ i ty an d m odel atmosphere s 
such as iHubenv fc Hubenvl ()1997f ) or iLi et al.l (|2005f) . 
iMortonson et al l (|2005r have argued that microlensing 
essentially measures the half-light radius {R1/2 = 2.44i?A 
for our model), so the first step should be to try to si- 
multaneously match our size estimates and the observed 
fluxes assuming this to be the case, since the many addi- 
tional parameters of the full disk models will make their 
direct inclusion in the microlensing calculations a major 
computational challenge. 

We thank E. Agol, M. Dietrich, C. Onken, B. Peterson, 
M. Pinsonneault, R. Pogge and P. Osmer for discussions 
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